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Analytic  Representation  of  the  Jacchla  1977  Model  Atmosphere. 

1.  INTRODUCTION 

The  diffusion  equations  for  the  Jacchia  1977  (J77)  density  model* 

cannot  be  integrated  in  closed  form.  Hence  numerical  results  must  be 

tabulated,  requiring  extensive  storage,  if  all  molecular  constituents  are 

tabulated  individually.  Fortunately  this  storage  requirement  can  be 

reduced  considerably  if  one  notes  that  the  solutions  for  constituents  in 

2 

diffusive  equilibrium  can  be  expressed  as  : 


•  exp  pi.  F  (hQ,  h,  T^j 


where 

ni  *  density  of*  i ^  molecular  constituent 
h  =  height 

T®  -  exospheric  temperature 
T  =  local  temperature 

ai  x  i*^1  constituent  thermal  diffision  coefficient 

M.j  =  i^  constituent  molecular  mass 

r  =  9(z)/[R*T(z,T«)]  dz 

0 

g(z)  =  acceleration  of  gravity  at  height  z. 

R*  =  8.31432X103  J (Kg-mol )'*/K 

In  general  one  then  need  tabulate  only  F  as  a  function  o.  h  and  T=o  and 
the  n.  at  the  boundary  hQ.  If  hQis  chosen  to  be  the  homopause, 
100km,  only  one  n.  need  be  stored,  from  which  the  others  are  derived. 

In  this  paper  we  explore  alternative  profiles  for  a  large  portion  of  the 
diffusive  region  (hM25km)  for  which  F  can  be  expressed  analytically, 
thus  further  reducing  storage  requirements. 
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2. 


TEMPERATURE  PROFILES  (h  >  125km) 


Since  a  fully  analytic  solution  for  the  temperature  profiles  specified 
in  J77  is  not  possible  it  seems  worthwhile  to  explore  other  forms.  In 
particular  the  Jacchia-Walker  (JW)  profile^: 

T lu(h,T  )  =  T  -(T  -T  )  e'aC 
leads  to  the  solution 

fl*.1.1,)'  9_ll 

R*  (Re  ♦  *0)2  T 

where 

E  =  (h  -  hQ)  (Rg  +  hQ)/(Re  +  h)  (aeonotential  altitude) 

Re  =  6356.766km 

9e  =  gravitational  acceleration  at  Earth's  surface 

=  9.80665  m/sec2 

T  =  T (h  ,  T  ) 

0  v  0  00  ' 

a  -  gradient  parameter 

For  application  to  J77,  it  is  prudent  to  choose  hQ  =  125km,  since 
this  is  the  inflection  point,  although  the  various  constituent  densities 
at  this  height  cannot  be  conveniently  derived  from  one  of  them,  as  at  100 
km.  The  gradient  parameter  may  be  chosen  to  fit  the  exact  J77  profile  in 
some  fashion,  such  as  to  match  the  slope  at  125km;  i.e. 

« =  V(T-  -  V  =  (dT/dh^T*  -  v 

Results  from  this  model  deviate  somewhat  fronr  thfe  exact  J77  modeTi  In 
particular,  for  an  exospheric  temperature  of  900K  the  local  temperature 
deviation  is  20K.  at  160km  and  the  mass  and  argon  densities  differ  by  6% 
and  12%  at  225km. 

To  obtain  better  fits  the  form 

T  >  [TJU  •  f]'1 


1 

r 

1  In’ 

T(h.  Tm  f 

oo 

-To 

> 

) 
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has  been  explored,  where  the  function  f  Is  chosen  to  retain  the 
closed-form  integrabillty  of  the  diffusion  equation.  In  addition  it  is 
desirable  that  '.he  resulting  total  temperature  fit  smoothly  the  exact  J77 
profile  at  h  -  hfl  =  125km  and  maintain  the  inflection  point  character 
(dT/dh  =  0).  It  is  evidently  the  lack  of  this  latter  characteristic 
which  makes  the  simple  JW  form  undesirable  and  possibly  accounts  for  much 
of  the  error  resulting  from  its  use.  These  modifications  lead  to  the 
following  boundary  conditions  at  £=  0. 

f  =  df/dS  =  0 
d2f/d  S2’  -%Z^Q2 

where 


G0W<,h]h  .  h 

0 

r,=  a  ♦  2/(Re  +  hQ)  =  G0/(T _  -  TQ)  ♦  2/(Re  ♦  hQ) 


For  convenience  the  functional  form  is  chosen  to  have  an  exponentially 
decreasing  character  with  parameters  hopefully  related  too  .  One  such 
form  is 


f  = 


0 

o 


C/(2B*Tj 


fS  M 
1 


where 


f8  (O  .  e'Sf>  (l  ■  e-^)r 

and  6  is  an  adjustable  parameter.  To  provide  further  flexibility  the 
form 


has  been  chosen.  The  additional  term  j  =  C„Cfp  satisfies 

2  *- 

2 

j  s  j'  =  j"  =  0 
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at  £  =  0. 


For  this  prof i le  the  solution  is  given  by 

f(h0.  h*  +  K{c-(  F  ^  q1  +  1/ 3 )/S 


JW 

+C 


>1  =  1  1  > 

pi  <V-at  -  bf>  s'  ‘ 


where  Fjw  =  Jacchia-Walker  solution 

Cl  ■  -c„  VIJB/V) 

a  =  -a ,  =  -1 ,  a  =-1/3 
i  J  ■> 

q,  =  e  Bit; 

bi  -  V’ 

K  »  9efte2/[R*(Re  ♦  hor] 
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1 

1 1/1 fij/B^ 2 ) 
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RESULTS 


The  FI etcher-Powel 1  non-linear  function  minimization  method  ’  has 

been  used  to  determine  the  parameters  B  ,  6  ,  c  to  minimize  the  sum  of 

1  2  2 


squared  residuals 


Q  ■  E(T  -  Tj?7)> 


for  T  =  500,  700,  900,  1100,  1300,  1500,  1700,  1900.  Bv  inspection  the 

co  * 

results  are  found  to  be  reasonably  well  represented  by: 


T.  <;  HOOK 

=  2o  1  +  10'“  (T  -  800) 


=  U.0215  -  0.005  ( T^  -  500)/200 

=  10~5l 0.0566  -  0.08X  +  0.04X2);  X  =  (T  -  900)/200 


1100K  <  T  <  1500K 


6  =  0.0385  -  0.012y  +  0.0l23y2;  y 


=  (T  -  1100)/ 400 


B2  =  0.0065  +  0.0167y 

C  =  10'5(0.0166  -  0.4548y7 ) 

2 


T  >  1500K 

Oo 

a  =  0.0388  -  0.0045z;  z  =  (T  -  15Q0)/400 

i  ® 

8  =  0.0232  -  0.0040z 

? 

C  =  -10'5  (0.4382  +  0.0387z) 


Mass  and  constituent  densities  were  computed  using  the  Jacchia-Bass  (JB) 
formulation  presented  here  for  the  T«  region  500K  -1900K  at  1Q0K  steps 
and  the  h  region  130km  -  lOCCkm  at  10km  steps.  The  results  were  compared 
with  those  from  the  exact  J77  temperature  profiles,  except  that  vertical 
flux  is  ignored  for  H.  Table  I  sutrcrtarizes  maximum  X  deviations  for  a 
6-constituent  gas  (O2,  0,  N^,  He,  Ar,  H)  with  the  total  mass  density 
given  by 

$ 

0  =  Z  H.  n./A;  A  =  Avoqadro's  Number 
i  - 1 

It  should  be  noted  that  a  true  J 7 7  calculation  would  call  for  different 
exospheric  temperatures  for  the  different  constituents  to  model  the 
different  diurnal  phases;  here  all  the  n^  are  computed  for  the  same 
T®  .  Argon,  with  the  largest  mass,  yields  the  largest  disagreement 
among  the  constituents.  The  agreement  is  quite  good,  well  within  typical 
model -experiment  differences  such  as  those  published  by  Forbes,  Marcos 
and  Gillette6,  and  Sharp  and  Brag.7  It  should  be  noted  that  the 
maximum  deviations  are  all  within  the  1100K  -  HOOK  region,  which  proved 
to  be  the  most  difficult  to  fit.  Outside  this  region  all  JB  mass 
densities  agree  with  J77  within  IX  at  all  heights.  Maximum  temperature 
deviations  are  13K  at  T„  =  1300K,  h  =  350km.  Ouside  the  1100K  -  HOOK 
region  the  maximum  disagreement  is  5K. 

A  FORTRAN  deck  encompassing  the  formulation  presented  here  is  available 
on  request. 


Table  I.  Maximum  Absolute  Value  of  0B-J77  Deviations  (X) 


TOTAL  DENSITY  Ar  DENSITY 

h  S  300 Km  1.1  (h  *  200,  =  1200)  1.5  (h  =  190,  T^  =  1 


h  i  500 Km  2.1  (h  *  500,  =  1300)  5.6  (h  =  500,  T^  =  1 

h  <  1000 Km  3.0  (h  =  750,  T^  =  1300)  8.5  (h  =  1000,  T^  = 


12- 


<1 


1 


4. 


CONCLUSIONS 


It  has  been  shown  that  analytically  solvable  temperature  profiles  can  be 
adjusted  to  realistically  represent  temperature  vs.  height  variations 
above  125km,  allowing  one  to  reduce  tabular  storage  requirements  and/or 
avoid  the  inconvenience  of  interpolation.  It  is  therefore  recommended 
that  similar  profiles  be  used  in  future  atmospheric  density  modelling. 
The  possibility  of  developing  analytically  integrable  profiles  below 
125km  will  be  studied  at  a  later  date. 
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